Electronic Structure of Ladder Cuprates 
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We study the electronic structure of the ladder compounds Sri^-zCa^Q^CUi and SrCu2C>3 . 
LDA calculations for both give similar Cu 3d-bands near the Fermi energy. The hopping parameters 
estimated by fitting LDA energy bands show a strong anisotropy between the t± and tu intra-ladder 
hopping and small inter-ladder hopping. A downfolding method shows that this anisotropy arises 
from the ladder structure. The conductivity perpendicular to the ladders is computed assuming 
incoherent tunneling giving a value close to experiment. 

PACS numbers: 71.10.Pm, 71.20.-b, 71.27. +a 



Sri4_ a; Ca 2: Cu2404iEla is the first material in which 
doped ladder (see Fig. |l|) can be experimentally stud- 
ied and compared to theoretical predictions of a 
Luther-Emery stat&, with a spin-gap, hole-pairing and 
superconductivityETE. A^pin-gap of A ~ 280K (for 
x = 9) has been measuredo and superconductivity under 
high pressure P > 3GPa has been founoHO in Ca-rich 
samples x ~ 1 1 (Call ) ,.. h aving ladder layers doped with 
20% holes (5 ~ 0.2)000. The transport properties are 
dominated by holes in the ladder planes. The normal 
state of the Call shows a strong anisotropy between the 
dc-resistivity along and .ac ross the ladder direction with 
p±/pn ~ 30 at T=100KE3. For lower temperature, both 
resistivities increase exponentially due to localization ef- 
fects. For higher temperatures, pn increases linearly in T 
while p± remains nearly constant, p± — 12mr2cm. The 
mean free path along the ladder is larger than the lattice 
constant while across the ladders it is smaller than the 
inter-ladder distance indicating incoherent transport in 
this direction. 

Moreover, fits of the spin susceptibility have shown a 
large difference between the exchaxiee coupling J±(J\\) 
along the rungs (legs) of the laddertSO even if both in- 
volve similar 180° Cu-O-Cu superexchange processes. 
Analysis of neutrpa scattering data gives J± — 72meV 
and Jy = 130meVlii For these reasons a detailed exam- 
ination of the electronic structure is desirable. 

In this paper we present LDA calculations of the elec- 
tronic structure which give estimates of effective hop- 
ping matrix elements between states on different copper 
ions. The LDA studies are performed for SrCu203, a 
compound which possesses the same kind of CU2O3 lad- 
der planes as Sri4_^GapCu2404i(see Fig. |l|). Recently, 
Arai and TsunetsuguEZI reported LDA calculations for 
M14CU24O41, (M=Sr or_Ca) which give similar results. 

The TB-LMTO ASAE3 energy bands for SrCu 2 3 are 
plotted in Fig. |^. The uppermost graph displays the 
bands on the path T = (0,0,0) Z' = (0,0,tt/2c), A' = 
(27r/a,0,7r/2c), and X = (27r/a,0,0). The two parallel 
bands near zero energy (Fermi energy of the half-filled 
band) are separated from the rest of the spectrum. They 
are due to hybridization through c-bonds of the 2p O- 



and 3d x 2_ z 2 Cu-orbitals. These bands hybridize with 
non-bonding O-bands near T. In the lower graph, the 
energy bands are shown up to the edge of the zone with 
Z = (0, 0,7r/c) and A = (2ir/a,Q,ir/c). Due to destruc- 
tive interference, the 2 D-bands do not display dispersion 
on the path ZA. 

The low energy physics can be described by an effec- 
tive model containing only these two bands with similar 
shape near the Fermi energy. Such a model includes only 
one state per Cu with effective hopping matrix elements. 
These 2 bands are the bonding (b) and anti-bonding (a) 
rung bands of the effective ladder model. 

Note that the parallel nature of the bands at k z = it /2c 
(path Z A') cannot be explained by effective interactions 
between nearest-neighbor (n.n.) Cu-sites only. In such a 
model, the dispersion along the k x direction is given by 
the inter-chain hopping matrix element between the sec- 
ond leg of one ladder and the first leg of the next ladder. 
Since the hopping matrix elements between two bond- 
ing(antibonding) states of the rung r of two neighboring 
ladders I and V b a ^ r {a a ,i, r ) = -^(4>a,i,r,i + {-)4><j,i,r,2) is 
given by 



(ba,l,rH(tlli)b ai l' tr ) 
{a>a,l,rH{tll')a a M ,r) 



7,Lr,2H(tll')(f>aM,r,l), 

^ CT ,i,r,2-H"(iH')^,'',ni)! (!) 



where the last index of (j> a ,i',r,i labels the ladder leg. 
Therefore b- and a-states should have an opposite disper- 
sion in the fc^-direction. Thus, an effective atomic model 
must contain some longer range inter-ladder hopping to 
account for their parallel nature. 

Arai and Tsunetsugu introduced a simpler rung pa- 
rameterization of the band structure fitting the b- and 
a-bands separately, allowing n.n. and n.n.n. hopping lead- 
ing to the formtO, 

e(k) = eo — 2/iim cos(fc z ) — 2/jn 2 cos(2fc z ) — 

1 3 
[ihx,i cos(-fc z ) + 4/i_l,2 cos(-fc*)] cos(k x ). (2) 

The values they obtained for 8^4^24041 are in good 
qualitative agreement with ours for SrCu2C>3 (see Ta- 
ble ph. Note that the signs of the interladder hopping 



1 



parameter h± i does not change between an a-band and 
b-band contrary to the expectations from Eq. |l|. All hop- 
ping parameters apart from /iy ;1 are higher order in the 
Cu-0 (t p d) and 0-0 overlaps t pp and thus much smaller. 

To gain more insight we introduce a single parameteri- 
zation of both bands in terms of intersite hopping param- 
eters shown in Fig. || The solution of this tight-binding 
model is 



e±(k) = € + €n(k z ) + e_L,i(fcz) cos(k x ) 



(3) 



±yje±,3( k *) 2 + e ±A k z) 2 + 2e -L,30 z )e_L,4(fc z ) cos(k x ) 
where 

e||(^ z ) = — 2tiM co8(k z ) — 2tii 2 cos(2fc z ), 

k 3k 
£±,i(k z ) = -4t_u3 cos(y ) - 4*_l,7 cos(-^), 

£±A k z) = *J-,1 + 2 *-L,4 COs(fc 2 ) + 2t_L, 6 C0s(2fc z ), 

fe 3fc 
ex,4(*«) = 2i±, 2 cos(y ) + 2t±, s cos(— -). (4) 

The parallel nature of the bands at k z = ir/2 is recov- 
ered if the term in the square root of Eq. |^ is indepen- 
dent of k x . This obtains if ej_ i 4(7r/2) ~ 0, leading to 
t±,2, t±,5 — or to fj_,2 — £_L,5- In these cases the disper- 
sion at k z — it/2 is due to the e± i term being a function 
of tj_,3 and tj_,7- It clearly shows that an intersite model 
must contain at least the third order hopping term £j_ 3. 

The coupled ladder system has a glide symmetry given 
by the product of a reflection through the c-axis (see 
Fig. [5j) and a translation of half a lattice constant along 
the ladder. When this operation is applied twice it is 
equivalent to a translation of one lattice constant along 
the ladder. This implies that the energy band at k x — 
and k x = it can be parameterized through one parame- 
ter. Actually e±(k x = tt, k z ) = e±(k x — 0, 2ir — k z ) as 
can be directly checked from Eq. ||. This allows one to 
use a single function containing the information about 
all hopping parameters. This symmetry also implies the 
lack of dispersion in the k x direction at k z = tt as dis- 
cussed above, as well as the square root form containing 
intra- and inter-ladder hopping terms. 

Introducing a = 1(— 1) for k x = 0(tt) Eq. || reduces to 
the simpler form 



e±,cr(fcz) = e +ey(fc z ) +crej_,i(fc z ) 
±(e_i_,3(M + fCj. ( 4(A z )), 



(5) 



representing the 4 energy bands of the double-ladder sys- 
tem. Rewriting the bands from the rung form (g) to the 
intersite form ^ gives the values of Table ||. They are 
consistent with each other and emphasize the dominance 
of the n.n. intra-ladder matrix elements f y and t± w.r.t. 
the others. Moreover, they show surprisingly that these 
two parameters describing n.n. Cu-Cu processes differ 
from each other by ~ 35%..— ..— . 

Recently, Andersen et alrBEl introduced a systematic 
downfolding scheme to obtain an effective single (or few) 



band model capable of reproducing the details of the 
LDA bands close to the Fermi level. Effective hopping pa- 
rameters are calculated by performing the Fourier trans- 
form of the downfolded Hamiltonian H (k) — ► H(R), for 
|R| less than a cut-off radius i?o- This has the advantage 
among others that it allows the origin of the parameters 
in the effective single or few band model to be traced. We 
have applied this scheme to the LDA bands for SrCu203 . 
The anisotropy fii 1 ^ t_L,i is best understood by starting 
with an effective model including the 3d x 2_ z 2 , 4s Cu- and 
2p x , 2p z O-orbitals (dsp-model). The dominant parame- 
ters (> O.leV) for in-plane hopping are given in the first 
rows of Table III. Here r labels the O on the rung of one 



ladder while 1 labels either a Cu or an O on the upper leg 
of the ladder. Vectors e x = (a/6, 0, 0) and e z = (0,0, c/2) 
gives the translation vector from one O(Cu) to the neigh- 
boring Cu(O) in the respective direction. The notation 
is such that £d(i)p x (+e :c ) denotes the hopping between a 
3d Cu-orbital at 1 on the upper leg to the neighboring 



2p x O-orbital at 1 + e^,. The on-site t t 



(l)d(i) 



hopping is 



non-zero as consequence of the downfolding of all other 
bands in the absence of local four-fold symmetry. 

The on-site energies of the rung-oxygen e(r) is slightly 
larger than that of the leg-oxygen e(l) due to the lo- 
cal environments. The distance between Cu and O is 
ri = 1.98A \\z and r r = 1.92A \\x implying a larger hop- 
ping along the x-direction according to t p d oc 1/r 4 and 
explaining the anisotropy for the tpd's listed in first row 
of Table |Tn] Moreover hopping processes involving s- 
orbitals are large with non-negligible t s d hoppings. They 
will contribute to the t p d hopping in second and higher or- 
ders through paths like d—s~p, d—s—s—p, etc. By down- 
folding the s-orbitals, these processes will strongly renor- 
malize the t pd such that *d{i)p,(+e^) > td{i)p w {-e x )- Lastly 
downfolding the p orbitals gives the effective (d)-model. 
Results are given in third column of Tabic [ll| They are 
consistent with our previous results. The anisotropy be- 
tween tj_,i and tiM can now be seen to be a consequence 
of different t p d hopping due to the downfolding of the s 
bands. 

The exchange interaction J between spins are usually 
difficult to estimate a priori. In perturbation theory, J 
scales as t^ d but given the large value of t p d relative to 
(ep ~ £d) the results are not reliable. Our analysis shows 
that the energy difference (e p — €d) does not differ much 
between rung and legs but t p d does. A ratio Ji i / J y j < I 
is expected consistently with previous resultsEj~E3. This 
will disfavor the hole coupling and the RVB liquid state 
and may explain why holes are unbound for T > 100K, 
however this point will not be discussed further. 

In the following, two limiting estimates of the con- 
ductivity are investigated. First, the band structure 
model limit is considered ignoring magnetic interactions 
although they must be important. In a metallic ground 
state, the conductivity reduces to an integral over the 
Fermi surface. 
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2e 2 
(2 



^ J dkr(e(k)) Vl (k)v 3 (k) (- 



df 

dc 



(6) 



Considering an electron filling of n = 0.8 the ratio be- 
tween the conductivity perpendicular and parallel to the 
ladder can then be simply computed yielding very large 
values i.e. <j\\/<t± — 75, 90, and 104, for the parameter in 
Table |j[ Thus roughly one has 



a»/a± ~ 100. 



(7) 



This large anisotropy is a consequence of a warped Fermi- 
surface with very small Fermi velocity perpendicular to 
the ladders. 

Second, as discussed previously, the resistivity data 
perpendicular to the ladder indicates a mean-free path 
smaller than the inter-ladder spacing. Holes seem to 
be confined to ladders and to hop incoherently between 
them. A detail structural refinement of CalltJ shows a 
complex distortion pattern of the Cu02-chains leading 
to substantial potential variations between rungs due to 
the proximity of apical O-ions at ~ 50% of the rungs 
An estimate of the interladder conductivity can be made 
by considering the limit where ladders form quasi one- 
dimensional metallic systems weakly coupled to each 
other. The conductivity is thus a consequence of the 
inter-ladder hopping term 



H = T k c\ k c 2 ,k + h.c. 



(8) 



where Tk = J2i hx,i(h)- The conductance a± is given by 

ot=A^N{er) 2 \T{k F )\\ (9) 

with M(ejr) denoting the density of state at Ep in the 
quasi-one dimensional metallic system. Results for the 
different parameter sets are a± ~ 0.1,0.09, and 0.074, 
respectively. They are thus all close the value of 



0.' 



mfi 1 cm , 



(10) 



which corresponds well with the experimental result at 
T > 100K. 

In this paper, estimates of the hopping matrix elements 
based on LDA calculations gave three main results. First, 
the effective intra-ladder hopping between n.n. are not 
the same, t\\ ^ t±. This was explained as the consequence 
of anisotropic t p d in the (pd)-model due to effective hop- 
ping through paths involving Cu s-states. Second, inter- 
ladder hopping is much smaller than intra-ladder and 
longer range hopping must be included. Third, estimates 
of the conductivity in the model where holes are unbound 
and confined into the ladder give good agreement with 
the experiment at temperatures T > 100K. 

One of us (T. M.) thanks the "Fond National Suisse" 
for financial support. We would like to thank O. K. An- 
dersen and M. Troyer for useful and fruitful discussions. 
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FIG. 2. LDA band calculations on the path TZ'A'XT (up- 





permost 


graph) and the path YZAXV (lower graph). 
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FIG. 1. The Cu 2 3 ladder plane. 



FIG. 3. Coupled ladders illustrating effective hoppings be- 
tween Cu sites (t) or Cu-Cu rungs (h) 
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eo 


-0.31 


0.46 


-0.44 


0.35 


h U 


0.41 
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0.45 


0.68 


ha 


0.08 


0.07 


0.08 


0.07 


h_L,i 


0.07 


0.03 


0.07 


0.03 


h±,2 


0.00 
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TABLE I. The hopping parameters in eV for 
Sri4Cu2404i and SrCu203 in the rung scheme. 



z' A' x r 




z a x r 





Sri4Cu2404i 


SrCu 2 3 


Fit 


Downfolding 




0.075 


-0.045 


-2.476, 




0.500 


0.565 


0.537, 


*X,1 


0.385 


0.395 


0.351, 


t±.2 


0.040 


0.040 


0.018, 


t±,3 


0.050 


0.050 


0.050, 


t±A 


-0.090 


-0.115 


-0.124, 


t-L,5 


0.040 


0.040 


0.053, 


*ll.2 


0.075 


0.075 


0.106, 


t±.6 


0.005 


0.005 


0.012, 


t±,7 


-0.020 


-0.020 


-0.023, 



TABLE II. The parameters in eV from direct fitting and 
downfolding of the LDA bands. 
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£rf e p ( r ) 


e p(l) td(\)pz(+ez) 


i(Z(l)p x (-e x ) 

\ /f X \ X 1 


^d(l)p x (+e x ) 


^Pi(l)p^(+e x +e,) 

f z \ ! f z \ t X i s / 


(dsp) 


-3.82 -3.66 


-4.00 0.74 


0.85 


-0.85 


0.30 


(dp) 


-3.83 -4.59 


-4.53 0.65 


0.51 


-0.63 


0.22 






^p z d')p,(+e,+e z ) 


(Dv-r- ( — e T +e, ) 

/ , VrS \ J- 1 Z / 


/■ -< V / 1: X \ "I I Z / 


^p z (l)p z (+2e.) 


(dsp) 


0.30 


0.01 


-0.13 


0.28 


-0.42 


(dp) 


0.62 


0.62 


-0.64 


0.19 


-0.05 




*Px(l)Px(+2e 2 ) 


i Px(r)p :c (+2e x ) 


e s 




i s (l)d(+2e z ) 


(dsp) 


-0.17 


-0.62 


2.43 


-0.24 


-0.14 


(dp) 


-0.09 


+0.31 












*s(l)d(-2e x ) 


*s(l)s(-2e x ) 




*s(l)p x (±e x ) 


(dsp) 


-0.32 


0.20 


-0.44 


1.89 


±2.10 



TABLE III. Hopping parameters in eV from the downfolding method. Notation is explained in the text. 
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